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ABSTRACT 

We present a theoretical study of CS line profiles in archetypal hot cores. We provide estimates 
of line fluxes from the CS(l-O) to the CS(15-14) transitions and present the temporal variation of 
these fluxes. We find that i) the CS(l-O) transition is a better tracer of the Envelope of the hot 
core whereas the higher- J CS lines trace the ultra-compact core; ii) the peak temperature of the 
CS transitions is a good indicator of the temperature inside the hot core; Hi) in the Envelope, 
the older the hot core the stronger the self-absorption of CS; iv) the fractional abundance of CS 
is highest in the innermost parts of the ultra-compact core, confirming the CS molecule as one 
of the best tracers of very dense gas. 



Subject headings: astrochemistry — line: 
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Introduction 
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Understanding how stars form at various red- 
shifts is crucial in order to infer how larger struc- 
tures such as galaxies are made and evolve in 
the Universe. To understand the process of star- 
formation, it is essential to determine the proper- 
ties of the gas in star-forming regions (hereafter 
called SFR), both in our own Galaxy and in ex- 
ternal galaxies. 

SFR encompass a large range of physical and 
chemical conditions. Within SFR, the gas and 
dust are recycled from prestellar cores to hot cori- 
nos and hot cores. During the star formation cy- 
cle, the pressure, density, temperature and chem- 
istry vary. 

Numerous papers present atomic and molec- 
ular observational surveys of prestellar cores, 
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Determining the gas 



and dust temperatures, gas density, molecular 
abundances, etc. of each component in SFR can 



only be achieved by a close comparison between 
observations and detailed modelling. 

In this paper, we present a theoretical study 
of the 12 C 32 S (hereafter CS) molecular emis- 
sion from hot c o res, m otivate d by the work of 
Dotv fc Neufeidl (|l997h and IMillar fc Hatchell 
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(|1998I ). Subsequent papers will present results 
for molecules such as methanol, HCO + and HCN. 
Our first aim is to provide observers with theoret- 
ical CS profiles that may help them interpreting 
the CS line e missions ar i sing f rom hot cores, such 
as those from Wu et al. I (l2010h . We do not specif- 
ically model any particular hot core. Instead, 
we model an archetypical hot core composed of 
ultra-compact core and surrounding envelope. We 
present estimates of line fluxes and line profiles 
for comparison with observations. 

Th e study performed by IMillar fc Hatchell 
(|1998l ) of the hot cores distinguished two zones of 
emission: the ultra-compact core (hereafter UCC) 
and the Envelope of the hot core. While the UCC 
zone is characterized by a size o f about 0.03 pc, an 
avera ge age of 3.2 xlO 3 yrs (see IMillar fc Hatchell 
19981 ) and an average density of lxlO 7 cm -3 , the 
Envelope corresponds to a more extended region 
(0.15 pc size) at a lower densi ty (lxl0 6 cm~ 3 , 
see again IMillar fc Hatchefllll998l ). In this paper, 
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studies of both the UCC and the Envelope CS line 
emissions are performed between lxlO 3 yrs and 
lxlO 5 yrs. It is expected that the Envelope sur- 
vives longer than the UCC once the protostar(s) 
is (are) formed. Thus, its emissions should remain 
detectable at later times than the emissions com- 
ing from the UCC. This is why, for the Envelope, 
we have investigated time up to 1 x 10 6 yrs. 

The key questions we aim to answer in this pa- 
per are i) what are the contributions of these two 
emitting zones to the total CS line emissions de- 
tected in hot cores? ii) How do these contribu- 
tions evolve with the age of the hot core? An- 
swering this question is crucial for improving our 
current understanding of massive star formation. 
CS is particularly useful for observers, as it emits 
quite strongly, not only in hot cores in our own 



Galaxy (e.g. i Beuther et al. 2002 ; Leurini et al 



see e.g. 



2007t lWu et alil2010ft b ut also in external galaxies 



Martin et all [2006; Bav et et al.ll2008a : 



Aladro et al 



20091 : iBavet et all l2009h . It is also 
recognized as to be one of the best tracers of 
very dense and warm gas with lin e critical den- 
sities of about 10 6 ~ 7 cm -3 dPlume et all Il992 : 



Linke k Goldsmith! 19801: ISnell et al.lll984h . In ad 



dition, its spectroscopic characteristics are very 
well knowrtj. 

In Sect. [21 we describe the models we use 
and present, for the first time, an interface code 
we have built between the UClH chemical model 
(hereafter called UCL_Chem) and the radiative 
transfer code SMMOL. In Sect. [3l we specify the 
parameters used for this particular study of the 
CS molecule emission in hot cores. In Sect. [4] 
we present our results and we show how the CS 
line profiles and fluxes vary with evolution of the 
hot core and different parameters. We discuss the 
results in Sect. [5] and conclude in Sect. [51 

2. Model descriptions 

To generate the CS line profiles, we have de- 
veloped an intuitive and friendly interface able to 
couple the UCL_Chem model and the radiative 
transfer code SMMOL. This Interface code will be 
very shortly publicly available. The UCL_Chem 
model is briefly described in Subsect. 1 2. II whereas 
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a summary of the main characteristics of SMMOL 
is presented in Subsect. 12.21 The interface is de- 
scribed in Subsect. 



2.1. The UCL_Chem model 



The UCL_Ch e m mo del is fully described in 
Viti k Williams! (1 19991) and its u pgrades are 
present ed in Viti et al. ( 20041 ) and IBavet et al 
(|2008bh . 



The UCL_Chem is a time-dependent 1-D chem- 
ical code that can be used to model the evolution 
of the gas and dust during the formation of a star. 
Here we use d it to simulate th e formation of a hot 
core. As in Viti et al. I (|2004h . we first model the 
collapse of a 10K core (Phase I); we then follow 
the chemical evolution of the remnant core once 
the star is born (Phase II). The presence of an 
infrared source in the centre or in the vicinity of 
the core is simulated by subjecting the core to an 
increase in the gas and dust temperature. 

In both phases, the chemical network is based 
on more than 1700 chemical reactions taken 



from the UMIST database (Millar et al. 1997 



Le Teuff et al"1l2000h involving about 200 species, 
of which 42 are surface species. The relevant sur- 
face reactions included in this model are assumed 
to be only hydrogenation, allowing chemical satu- 
ration where this is possible. 

One of the outputs of the UCL_Chem is the 
fractional abundance (with respect to the total 
number of hydrogen nuclei) of gas and surface 
species. See Sect. [3] for a description of the grid 
of UCL_Chem models ran for this study. 

2.2. The SMMOL model 

The molecular line radiative transfer code 
solves the multilevel radiative transfer problem 
in a 1-d spherical geometry. The code we use, 
SMMOL, is base d upon two codes; Multi-Mol 
(lYates et al.lll997l) an d the SMULTI code devel- 
oped bv lHarperl (|l994h . SMMOL uses an Acceler- 
ated Lambda Iteration (ALI) scheme to speed the 
convergence of the iterative scheme that is used 
to solve a set of linearly perturbed kinetic master 
equations, in order to determine the steady state 
populations of a molecule's energy levels and the 
radiation f ield. This is the MULTI m ethod de- 
scribed in IScharmer k Carlsso n (1985). Subse- 



quently the MALI method ([Hummer k Rvbicki 
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19921 ) was added to SMMOL; this uses an ALI 
technique to speed the convergence of a set of 
preconditioned kinetic master equations. 

SMMOL (Spherical Multi-MOL) is a gen- 
eral non-LTE molecular line radiative trans- 
fer code that has reproduced the spectral lines 
observed towards a large number of sources 



e.g. iBenedettini et ail 120061: lYao et all 120061: 



Lerate et al.l2008T Tsamis et al.||2008l;lLerate et al 



2010). It is fully described in iRawlings fc Yates 
(|200lh and has been bench marked with other 
radiat ive transfer codes in Ivan Zadelhoff et al. 
(|2002l ). We recommend these papers to any 
reader. The output line profiles are convolved 
with user supplied t elescope properties using the 
method described in Schoenberg ( 19881) . 

Typically SMMOL uses 400 rays to compute 
the spherical cloud modelled to compute the inten- 
sities at each radial grid-point (hereafter " shells" ) 
in the cloud. The code is capable of adapting its 
sampling along each ray to take account of large 
velocity changes between shells e.g. if the line 
of sight velocity change between adjacent shells 
causes the individual line absorption profiles at 
these adjacent shells to be non overlapping in fre- 
quency space; this is not physical and can allow 
photons to escape the cloud that would otherwise 
have been absorbed. 

2.3. The Interface code 

The Interface is a fortran 95 programme, that 
transforms the output from UCL_Chem into line 
fluxes and profiles via the use of SMMOL. It is 
able to model various gas phases from diffuse gas 
to hot cores. It is currently being developed for 
AGB stars and planetary nebulae. 

The output file of the UCL_Chem model (i.e. 
fractional abundances as a function of time and 
depth) is the input file of the Interface. The 
UCL.Chem provides a grid in optical depth (Av) 
of fractional abundances of about 200 species, at 
various time steps. At each time step, the grid in 
Av has to be adapted to the spatial (linear) grid 
of shells used later on in SMMOL and described 
by: 

Av x (1.6 x 10 21 /met) .. .„ „ ,. 

n = H ; '- '- - (Ar/2.0 + (Ar x i)) 

density 

(1) 

where met is the metallicity assumed in the 



UCL_Chem and Ar is the thickness of the shells 
(assumed all equal here - linear distribution of the 
shells) : 

size , „ 

(2) 



Ar 



n s hell 



Then, the Interface also manages the allocation, 
to each shell of the appropriate physical values, 
required to run SMMOL. They are the density, 
temperature, fractional abundance etc which are 
extracted from UCL_Chem. To do so, the Inter- 
face converts Av into distance and interpolate the 
UCL_Chem values alongside the shell grid using 
SPLINE and SPLINT function^. For this study, 
the dust temperature is assumed to be equal to the 
gas temperature since in the UCL_Chem model, 
the dust temperature is not calculated. For our 
models (see Sect. [3]), this assumption is valid since 
the opacity and the density are very high in both 
the Ultra-compact core and the Envelope (Av 3> 
100 mag) . Once the Interface has created the cor- 
rect grid, it automatically runs the SMMOL pro- 
gramme. The outputs of SMMOL are plotted and 
tabulated by the Interface. 

We hope that the Interface will be used to 
automatically interpret data in from space tele- 
scopes (e.g. Herschel, JWST) and observatories 
(e.g. ALMA, e-Merlin...). 

3. Choice of parameters 
3.1. General assumptions 

The UCL_Chem models are converted into in- 
puts for the SMMOL code using the Interface 
code, as described in Sect. 12.31 We describe here 
our choice of parameters and assumptions. 

The main sampling parameters such as the 
number of radial density shells, and the line-of- 
sight frequency sampling through the cloud, are 
determined beforehand to ensure that the output 
fluxes are invariant with respect to sampling for all 
the models run through SMMOL. We found that 
n s hell — 100 and 400 lines-of-sight for ray trac- 
ing, ensured that the results were invariant; the 
smallest clouds we sample are actually very over- 
sampled. 

We used the first 40 rotational levels of CS in 
the vibrational ground state; the molecular data 
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and the collisional rates with respect to to H2 are 
from the Cologne Database for Molecular Spec- 
troscopjQ. 

The kinetic temperature law (see Eq. [3]) we 
used is a compromise between the need for a 
power law which is flatter in the inner hot core 
(to take into account photon trapping effects due 
to higher optical depths, which slow the cool- 
ing of material by radiation) and the need for a 
1/r 5 law to describe the cooling we would see 
in the lower optical depth outer cloud. We have 
use d Eq. [3] to be consistent with previous stud- 



ies dMillar fc Hatchelllll998l:IViti fc Williamsll 19991: 



Benedettini et al.ll200d : iLerate et al.ll2008ll201(£ 
We chose a turbulence velocity of 1.5 kms -1 
(since observations sh ow narrow line profile as in 



Hatchell et al.lll998bl ) assumed no velocity gradi- 



ent, a typical distance of 450pc (Orion-KL: see 



de Vicente et al.l |2002) to the source and solar 



metallicity. If these hot cores had a large velocity 
gradient, it would reduce the optical depth of the 
lines and allow photons to travel more freely in the 
hot core. The spectra will therefore become less 
absorbed and line broadening due to optical thick- 
ness would be reduced. However the line width 
would be increased as emission would come from 
a greater range of velocities. Actually, there is cur- 
rently little evidence for velocity gradients in these 
systems. Observations of lines from these sys tems 
show narrow lines (e.g. lHatchell et al.|[l99 8b') sug- 
gesting that all parts of the UCC and Envelope can 
be radiatively coupled. The line absorption pro- 
file is clearly dominated by turbulence, with the 
kinetic temperature contributing to the FWHM. 

The cloud was illumina ted by the sta ndard in- 
terstellar radiation field ( Habinel 19681) an d the 
dust extinction model is from Mathis ( 1990) • The 
emergent spectra were convolved with the appro- 
priate telescope beam (either JCMT or IRAM). To 
enable us to make useful predictions for ALMA or 
any interfcrometric observations (e.g. CARMA, 
PdBI,...), we need to predict emission from the 
source at sub-arcsec resolution. The 24 models 
whose parameters are displayed in Table Q] have 
been run twice, once with IRAM/ JCMT resolu- 
tions (line profiles presented as such, see Fig. [2J- 
2]), and once with effectively an infinite resolution 
(radial distribution, see Fig. [5]) allowing us thus 



to get the largest range of (sub-arcsec) resolutions 
possible. The gas density (see Eq. |4|), fractional 
abundances and cloud size are provided by the 
UCL.Chem model (see Fig. [[]). 

Currently the high resolution observational 
data that exist suggest that objects are ellip- 
soidal or spherica l in shape ( Davis et al. 2010; 
Graves et allbOlOl IWu et alj|2010h . We are con- 



strained by the 1-D nature of the UCL_Chem and 
SMMOL models, which means the only models 
we can construct are spherical. The most likely 
effect of non-spherical geometry would be the sce- 
nario where collapse was aided by magnetic field 
lines, giving a flattened density profile with an 
equatorial enhancement of material. The incli- 
nation of the source now becomes important; a 
pole-on source will have lower optical depth than 
an equatorial-on source. This means we could, for 
instance, overestimate the number of low optical 
systems, and underestimate the mass of objects. 

We ran the Interface code for several ages 
from 1 x 10 3 yrs to 1 x 10 6 yrs (see Table [T|). 
iMillar fc Hatchelll (|l998h assumed more specifi- 
cally a typical age of 3.2 x 10 3 yrs for the Ultra- 
Compact Core (see Models HCi, HC3, HC5 and 
HC7 to HC13) whereas a longer time for a less 
dense gas such as the one contained in the En- 
velope (see Models HC2, HC4, HC6 and HC14 to 
HC24) is expected. 

3.2. UCC and Envelope specific assump- 
tions 

To reproduce the CS line emission in hot cores 
for various density and temperature structures and 
ages, we have run over 80 UCL.Chem modelt@but 
present here only results from the most interest- 
ing ones; the parameter choices made for these 24 
calculations are summarized in Table [TJ 

The size of the region has been set to 0.03 pc for 
the Ultra- Compact Core and 0.15 pc for the En- 
velope. Other input parameters such as the FUV 
radiation field, the cosmic ray ionisation rate, the 
gas-to-dust mass ratio,... are all set to their stan- 



*See the CDMS website: |http: / /www. astro. univ-koeln.de/cdms/] 



3 The grid of models can be found at 
http:/ /www. homepages. ucl.ac.uk/^ucapdwi/interface/ 
Still under construction, this website aims at providing 
a variety of UCL_Chem models that the user can couple 
with the radiative transfer code SMMOL to obtain, for a 
various range of physical and chemical conditions, the line 
intensities and line profiles of more than 200 species. 
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dard values as in lBavet et al. | (|2Q08bl ) (see Tables 
1, 2 and 3 in their paper). 

To represent the UCC zone, we ran the 
UCL_Chem model with a temperature of 10 K 
in Phase I, collapsing the core to a critical density 
of 10 7 cm -3 (central density: n c ). We chose a cen- 
tral density of 10 7 cm -3 in order to b e consi stent 
with our previous work iBavet et all (|2008tJ l and 
because it is the density derived by single-dish ob- 
servations (H atchell et al. 1998a). In fact, the cen- 
tral density of hot cores ma y be higher than tha t 
(see interferometric data in iBeltran et al.l [2005) . 
However due to the high opacity of hot cores (Av 
;2>100 mags) a small change in the central density 
should have a negligible effect on the fluxes. 

In both regimes (i.e. UCC and Envelope), when 
no temperature and density profiles are applied 
(see Models HCi and HC2), the temperature and 
the density are kept fixed at 300 K and 1.0 x 10 7 
cm -3 , respectively, for the UCC zone, and to 101.5 
K and 3.8 x 10 6 cm -3 , respectively, for the En- 
velope case. For Models HC3 and HC4, we only 
applied a density profile (see Sect. 14. 1|) . 

When a temperature profile is applied (this 
was done for Models HC5 and HC7 to HC13 us- 
ing Eq. [3]), the temperature varies from 120 
K to 300 K (between the lowest to the highest 
Av, respectively ) . We used the formula seen in 
Viti fc William) (|l999h : 



T(r)=T c x(-)-° A 
ro 



(3) 



where T C =300K is the central temperature typ- 
ical for hot cores (see iMillar fe Hatchelll Il998l 



Viti et al.l 120041 ) and r =0.18 pc is the distance 



from the edge of the hot core to the newly born 
star (see Fig. [1]). When a density profile is applied 
(e.g. Models HC 3 , HC 5 and HC 7 to HC13 and Eq. 
|4|) the density profile, from the lowest to the high- 
est Av, leads to the density varying from 4.9 x 10 6 
cm~ 3 to 1.0 x 10 7 cm" 3 respectively. We use d 



the same formalism as in ICaselli fc Myers] (Il995|) ; 
lHatchell et"ail (|2000h and lBacmann et alj (|2000), 
which is effectively a Bonnor-Ebert sphere approx- 
imation: 



n(r) = n c x (H ) 



-1.5 



(4) 



(jNomura fc Millarl 120041) is the radius between 
isothermal and non-thermal velocity effects in hot 
cores. 

To model the Envelope, we have run the 
UCL_Chem model with a temperature of 10 K 
in Phase I, letting the UCC collapse to a critical 
density of 3.8 x 10 6 cm -3 (i.e. the density of the 
inner shell of the Envelope, see Fig. [1]). From 
the lowest to the highest Av, when a temperature 
profile is applied (e.g. Models HCg and HC14 to 
HC24 using Eq. [3]), the temperature varies from 
58.5 K to 101.5 K, respectively. When a den- 
sity profile is applied (e.g. Models HC4, HCg and 
HC14 to HC24 using Eq. |4]) the density varies from 
1.0 x 10 6 cm -3 to 3.8 x 10 6 cm -3 , respectively. 

4. Results 

In the following section, we present the results 
for our study to note the effects of the variations in 
the internal structure of the hot core (see Subsect. 
I4.ip . and its age (see Subsect. 14. 2p . and we com- 
pare the spectra from the UCC and the Envelope 
(see Subsect. 14.31) . having addressed velocity and 
geometry effects in Sect. 13.11 Figures [50 show 
examples of the most interesting changes affecting 
the line profiles whereas Table Q] summarizes for 
each hot core model the integrated line fluxes of 
the CS transitions obtained. 

4.1. Influence of the density and tempera- 
ture profiles on the CS line emissions 

The influence of the density and temperature 
profiles on the CS line emissions is shown in Figf2] 
(from bottom to top). This figure represents thus 
results obtained for Models HCi, HC3 and HC5 
(for the UCC) and Models HC 2 , HC 4 and HC 6 
(for the Envelope) whose parameters are seen in 
Table [TJ 

In the case of constant density (top plots of 
Fig. [2]), it is interesting to note the differences 
in the line profile shapes of the low-J CS lines (up 
to CS(3-2)) as compared to those of the higher- J 
CS transitions (from CS(4-3)). Firstly, high-J CS 
lines show the strongest emissions. Secondly, the 
low-J CS lines have a narrower line width than the 
high-J CS transitions (by a factor of about < 1.5- 
2.0). The first result can be understood by look- 



where n c — 10 7 cm 3 is the density assumed at 
the center of the hot core and r c = 0.05 pc 
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derived from the models 
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ing at the level population distribution. Indeed, in 
Model HCi, the majority of the collisions occur in 
the high levels of CS, favoring transitions at high- 
J rather than low- J, whose levels are less popu- 
lated by one order of magnitude on average. In 
addition, the transitions between the high-J levels 
have high A^ coefficients {A^ is proportional to 
is 3 ). These factors give also rise to higher source 
function terms and so to brighter emission. In 
parallel, the high-J lines, as well as being brighter 
than low-J transitions, are also broader and show 
more spectral structure than the low-J transitions. 
Indeed, the higher-J lines have flatter peaks and 
some show strong self absorption at the systemic 
velocity. The line broadening is a consequence of 
opacity, more often seen in stellar spectra and of- 
ten called the curve-of-growth. The low optical 
depth lines have narrow Gaussian line shapes. As 
the optical depth increases the line centre emis- 
sion saturates, however the line wing emission can 
still increase and it is this increase that broad- 
ens the line and increases the FWHM of the line. 
Eventually the line wings saturate and flattened 
line-shapes are observed. Finally the optical depth 
at line centre is high enough to promote self ab- 
sorption and twin-peaked spectra are formed. The 
observed molecular line curve-of-growth spectral 
behavior can give a spectral signature for warm 
dense hot cores. 

To disentangle the influence of the density and 
the temperature profiles on the CS line emissions, 
we have first kept the temperature constant to 300 
K in the UCC and 101.5 K in the Envelope what- 
ever the UCC (Envelope) radial shell (see Models 
HC3 and HC4, respectively). In such case (den- 
sity profile only), the UCC CS line profile shapes 
do not seem significantly changed as compared to 
the case where there is no profile. On the contrary, 
for the Envelope, where the difference in densities 
between the outer and the inner shells is steeper 
than for the UCC, we see more significant changes. 

However, one notes that for Model HC3, a slight 
broadening of all the lines is seen, between a factor 
of 0.33 and 0.16 for the CS(l-O) and the CS(7- 
6) transitions, respectively. In addition, the peak 
antennae temperatures are on average weaker than 
in the case where no profile is applied (differences 
varying between 3 K and 12 K). A saturation in 
the CS(7-6) profile is also seen. This is due to the 
higher fractional abundances of CS obtained from 



the chemical model (see Table [lj. 

Finally when we add the temperature structure 
(see Model HC 5 for the UCC and Model HC 6 for 
the Envelope), we see that (Fig. 2), the peak an- 
tennae temperatures of all the CS lines but CS(1- 
0) in the Envelope are indeed weaker by a factor of 
1.5. To be more precise, we found that the mod- 
elled Tpeaks change in all 80 models when we im- 



plement a 



temperature profile as compared 



to their values without. It means that the distri- 
bution of the temperatures inside the hot core, i.e. 
the temperature variation seen from shell to shell, 
do have a consequence for the integrated (" total" ) 
profiles of CS lines. The differences in the Tpeak 
values range from 20% to a factor of two, depend- 
ing on the line and the source considered (see in 
Fig. 2 the bottom two plots for UCC and Enve- 
lope: both show a decrease of the CS line Tpeak 
when a temperature profile is implemented). We 
believe that ALMA may be sensitive to factors 
as small as this in Tpeak. In fact already with 
the CARMA and IRAM-Plateau de Bure Interfer- 
ometer such fac tors are detectable for hot cores 
dWu et al.ll2oToh . It may therefore be possible to 
estimate Tc and potentially reconstruct the tem- 
perature profile of the observed source by using 
Eq. 3. 

All the lines (except J=l-0) are thus good 
probes of these changes in temperature. The peak 
antennae temperature is a measure of the gas ki- 
netic temperature if the gas is in LTE, the column 
is optically thick and the source is resolved by the 
telescope at the observing frequency. 

The CS(l-O) does not seem affected at early 
times by the changes in temperature in the En- 
velope. This might come from the fact that the 
fractional abundance of CS is quite low in the En- 
velope. The is also 300 times less for the J=l-0 
transition compared to the J=7-6 transition. Also 
the energy of the J=l level -E(j = i) = 2.35 K, which 
means that in a gas of 50-100 K its population 
is comparatively less than in clouds where T=10 
K. This is why the 1-0 line is the weakest line 
in both the Envelope and the UCC. However the 
fractional abundance of CS increases by 20 in the 
Envelope during its evolution and that is why the 
line strength grows and the line starts to become 
flat topped. 
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4.2. Influence of the age of the hot core on 
the CS line emissions 

Figure [3] shows the evolution of UCC CS line- 
shapes with the evolution of the hot core. We 
see the high-J lines increase in flux with time and 
become increasingly flattened, broader and self 
absorbed; some eventually produce twin peaked 
spectra. The low-J lines grow to large intensities 
with time and also show flattened profiles at large 
times. 

Figure U shows the evolution of Envelope CS 
lineshapes with the evolution of the hot core. It 
shows that transitions for J=5, 4, 3, 2 are initially 
the brightest transitions, with the J=2-l transi- 
tion being the brightest line; the J=3 to 7 transi- 
tions are self absorbed at all times, with the J=7-6 
line being the broadest line. The J=l-0 transition 
is initially the least bright and narrowest line. 

By the largest time, the J=7 to 2 transitions 
are now all self absorbed and have basically the 
same width. There has been a modest increase in 
the peak flux of these lines. The J=I-0 line is now 
very bright, broader and has a fiat top. 

The UCC lines are brighter than the Envelope 
lines. For the Envelope, the best tracer of age 
seems to be the CS(I-O) transition which shows 
the largest variations in fluxes with respect to 
time (see Models HC14 to HC24). In principle, 
this transition could be used as an evolutionary 
indicator. However we note that the sulphur- 
bearing chemistry depends very critically on the 
gas temperature and density (as als o found by 
IWakelam et"aDl2004t IViti et alJIiool and hence, 
care must be taken to constrain initial conditions. 
In the case of the UCC, the differences in fluxes 
from 1 x 10 3 yrs to 1 x 10 6 yrs may not be large 
enough to be detectable. 

4.3. Influence of the size of the hot core 
on the CS line emissions 

There are, as expected, clear differences be- 
tween UCC CS line emissions and those coming 
from the Envelope (see Figs. [3J [3] and [U and Table 
[T|). We restrict our comparison to the most real- 
istic cases which are the cases where both density 
and temperature profiles are applied to the two 
zones. 

A first interesting remark concerns the radial 
distribution of the CS fractional abundance, for a 



given time. From column 9 of Table [TJ for a fixed 
time (e.g. 3.2 x 10 3 yrs ) one notes that, from the 
inner UCC shell (closest to the star) to the outer 
Envelope shell (edge of the hot core) , the CS frac- 
tional abundances does not increase nor decrease 
constantly (see paired models in Table [TJ. In fact, 
Tabled] shows that the CS fractional abundance is 
predominantly produced in the UCC inner zone of 
the hot core, confirming CS as preferentially linked 
with very dense gas, and therefore could be consid- 
ered as an ideal tracer of this gas phase (se e also 
Linke fc Goldsmith! Il98lj ISnell et al.lll984l ). For 



example, we can see that for the coupled models 
HC 3 and HC 4 , the UCC (HC 3 ) abundance of CS 
is 1.72 x 10~ 8 whereas the Envelope (HC4) can 
account for only 3.77 x 10~ 9 . This is a factor of 
4.6 difference. In other words, the UCC produces 
4.6 times more CS (in abundance) than the Enve- 
lope. For other models, the difference are of 2.6 
(HC 8 and HC15), 7.2 (HC 9 and HCi 6 ), 9.6 (HC10 
and HC17), 12.1 HCn and HCi 8 ), 17.6 (HC12 and 
HC 20 ) and 11.5 (HC13 and HC 22 ). The only mod- 
els where the CS is not predominantly produced 
by the UCC is the pair HC7 and HC14 (difference 
of 0.9 only). 

The radial integrated CS line fluxes distribution 
is shown in Fig. [5] at 3.2 x 10 3 yrs (as assumed 
in lMillar fc Hatchelll[l998l ). We note that the line 
flux distribution as a function of radius at different 
times (not shown) are the same as those shown in 
Fig. [5] Note that, for the plots shown in Fig. [SJ we 
did not apply any telescope convolution when we 
have run SMMOL because we wanted to identify 
all the emission from the UCC and the Envelope 
(hence the differences in CS line fluxes between 
Table [Band Fig. [5]). For both plots in Fig. [SJ 
the protostar is located at the origin of the x-axis 
(right hand side) following the convention adopted 
in Fig[TJ The maxima integrated CS line fluxes are 
located on the plots by a thick black cross. The 
maxima are all distributed in the UCC within a 
restricted radius (8 — 9 x 10 16 cm from the source) 
whereas for the Envelope, their position is more 
spread (3.5 — 5.5 x 10 17 cm from the source). In 
the Envelope, the CS lines fluxes maxima, from 1- 
to 10-9 are distributed deeper and deeper inside 
the gas. A turnover occurs for the maximum of the 
CS(10-9) line flux. From this transition onwards 
up to the CS(15-14) line, the maxima location is 
moving towards the edge of the Envelope. The 
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turnover is actually controlled by the competition 
between the source function and absorption terms 
at each radius. These terms depend upon volume 
of gas at a radius, the population of the levels and 
the Einstein A coefficients of the transitions (Note 
that Aij increases with J). In the Envelope this 
produces most flux for the J=10-9 transition. The 
turnover at large radii is caused by the reduction 
in the source function and optical depth due to 
decrease of gas density and temperature. 

In the UCC, the same factors affect the line 
emission as for the Envelope. The higher den- 
sity and temperature explain why the J=13-12 
emission is the brightest transition. The turnover 
happens proportionately much closer to the UCC 
outer edge than in the Envelope and thus does 
not appear clearly. This is because the optical 
depth and source functions of the optically thick 
lines drops due to falling densities and tempera- 
tures and these actually cause the low J lines to 
turnover as well. 

Here, we have simplistically assumed that when 
a hot core is observed the total CS flux (taking 
into account the self-absorption) is, to a hrst or- 
der approximation, equal to the sum of the emis- 
sion coming from the UC C and from the Envelope. 
Millar fc Hatchelll (jl998l) in their study, assume 
that hot cores can be as large as 1 pc size and 
consequently that there is a potential third con- 
tribution toJhe_Jotal^nokjcuh2xemission from the 
halo fsee lMillar fc Hatchell|[l998l ). We did not in- 
vestigate such region in this paper since the halo 
is much more diffuse than the Envelope and the 
UCC, and the CS may not be a good tracer of 
such gas. 

As seen from Fig. [SJ contrarily to the CS(1- 
0), (2-1) and (3-2) line fluxes, it is interesting to 
note that the high-J CS line fluxes (from 6-5) are 
mainly emitted from the UCC zone (factors of dif- 
ferences between the Envelope and the UCC con- 
tributions from 3.6 to 13.1 - see also Table [1]). This 
makes the low-J and the high-J CS lines better 
tracers of the Envelope and the UCC, respectively. 
This result can only be confirmed by interferomet- 
ric observations. The envelope has a lower hydro- 
gen number density and kinetic temperature, i.e. 



E u= 



(J=7) 



49.3 K than the UCC. The CS rotational 



the excited column at high-J levels in the UCC is 
brighter than that in the Envelope. In the UCC 
(T=100-300K) the J=l level is very underpopu- 
lated as compared to the Envelope and therefore 
has a much lower optical depth and source func- 
tion; it is for this reason that despite the fractional 
abundance of CS increasing we see a weak J=l-0 
line at all times in the UCC. 

Finally we note that in the case of the Envelope, 
the lines between J=2-l and J=6-5 show similar 
widths (within a factor of 1-2 kms -1 ), contrarily 
to the UCC case where line widths are quite differ- 
ent from a transition to another. The high-J lines 
remain optically thin for much longer in the UCC 
than for the Envelope. The absorption profile of 
each line is slightly broader because of the higher 
kinetic temperature in the UCC (0.2-0.3 kms -1 ) 

5. Discussion 

Unfortunately, there are no complete (i.e. from 
J=2-l to J=7-6) datasets of interferometric CS 
observations for hot co r es. CS data presented i n 
Hauschildt et all (ll995l):lMacDonald et al.l (Il996h : 



Beuther et al l (200z 
Leurini et al.l ( 2007 




Beltran et al.l ((20051 ) and 
focuss only on certain tran- 
sitions of CS or on isotopologues of CS. On the 
other hand, multiples lines of 12 C 32 S are observed 
in galactic hot cores but only using single-dish 
instruments. Despite these limitations, qualita- 
tive com parison of our result s with the data pre- 
sented in lMurata et al.l (ll994T); IChandler fc Wood 
( 1997 ) and more recently in Wu et al. ( 2010l) can 



be made. Of particular interest, IWu et all (|2010h 



high-J levels in the UCC are therefore more highly 
populated that the same levels in the envelope. 
Although the UCC is smaller than the Envelope 



found that the CS(2-1) transition is indeed less 
compact than the CS(7-6) as seen in maps of 50 
massive very dense galactic sources. This agrees 
very well with our model predictions and seems 
to confirm that the high-J CS lines are one of the 
best tracers of very dense compact gas (i.e. better 
than HCN). Similarly, they found that their mean 
and median linewidth increase for high-J CS lines, 
which is also in agreement with our model predic- 
tions. Here we do not attempt to model any of the 
sources of their sample but the fact that our the- 
oretical approach and their observational results 
converge is encouraging. 

Fig. [5] thus used by observers to estimate in- 
tegration times as it gives the expected flux of 
the first 15 transitions of CS at various resolution 
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(i.e. various radii). The fluxes are generally pro- 
duced by a combination of a chemical model, line 
and continuum radiative transfer code and tele- 
scope convolution algorithm^]- The deconvolved 
fluxes presented in Fig. [5] can be used to estimate 
fluxes for other telescope parameters and source 
distances. For different source types and masses, 
an observer can scale these results, but the scaled 
results would be approximate at best. With in- 
terferometers such as the IRAM-Plateau de Bure, 
spatial resolutions up to 7 x 10 16 cm are already 
accessi ble (for a distance of 4 50 pc). As an ex- 
ample, [de_Vlcente_£t_alJ (|2002l ) already performed 
a detailed study of HC3N in the Orion KL hot 
core. With Fig. [5J we show that similar studies 
are possible for CS. Indeed, a reso lution of 5 " 
is reached in lde Vicente et al.l (|2002l ) work, which 
corresponds to 0.01 pc (i.e. 3.856 xlO 16 cm) i.e. 
the CS UCC zone emission. More information on 
the structure of the CS emission in hot cores will 
be obtainable soon with ALMA. 

6. Conclusions 

We have performed a systematic theoretical 
study of the line intensities and the properties 
of CS in archetypical hot core environments. We 
have coupled via a user-friendly interface, a large 
grid of chemical models with SMMOL radiative 
transfer code and obtained line profiles of the 
CS(l-O) to CS(15-14) lines for a variety of den- 
sity, temperature, size and age. We also provide 
(Fig. [5]) observers with estimates of line fluxes at 
various resolution for the first 15 transitions of CS 
molecule. 

Our main conclusions are: 

• the CS fractional abundance is highest in 
the innermost parts of the UCC whatever 
the age of the hot core. This confirms the 
CS molecule as one of the best tracers of the 
very dense gas component (see Sect. 14. 3[) . 

• the high-J CS lines have the strongest line 
fluxes, and the linewidths are broader than 
those of low- J CS lines (see Sect. 14. ip . 



The peak antennae temperature of all the 
CS transitions except for the CS(l-O) line is 
a good tracer of the kinetic temperature in- 
side the hot core because it is very sensitive 
to its changes (see Sect. 14. ip . 



• In the Envelope, the older the hot core, the 
stronger the self-absorption of CS. The best 
tracer of age seems to be the CS(l-O) line 
which show the largest variations in fluxes 
with respect to the time (see Sect. 14.21) . 



• The CS(l-O) flux is coming mainly from the 
Envelope while the high-J CS line fluxes are 
better tracers of the UCC zone (see Sect. 



7 We will publish on the web the deconvolved 
fluxes as well as the convolved fluxes. See 
http: / /www.homepages. ucl.ac.uk/~ucapdwi/interface/ 
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AT A TIME t : 



Envelope -0.15pc size (100 shells, 



UCC = 0.03 pc size (100 shells 




size 



z 100 = outer UCC shell 
n 100 =4.9 x10 6 cm" 3 
T 1oo =120K 
Y(CS) t> 

100 



z-,= inner UCC shell 
n-!= 1x10 7 cm" 3 
T-,= 300K 
Y(CS) t>1 



16 cm- 3 



z 100 - outer Envelope shell 
n 1oo =1.0x10 e 
T-/QQ- 58.5 K 

Y(CS)t,ioo 



i6 cm -3 



z 7 = inner Envelope shell 
n^Z.SxICf 
T 7 = 101.5 K 

Y(CS) ti1 



Fig. 1. — Cartoon representing the structure of the hot core we adopt: the Ultra-compact core UCC (right 
hand side) and the Envelope (left hand side). The characteristics of each zone, in the multi-points approach, 
for a certain time t, are specified. The grey scale does not have any physical meaning but should be considered 
more as an guideline for the eyes. It symbolises the gradient in density and temperature through the hot 
core with the darkest zones corresponding to the highest values of density and temperature. The shells of 
width Ar are represented by vertical bars from number 1 to 100 (zi showing the location of the shell number 
1). Regular and italic fonts are used for distinguishing the UCC shells from the shells of the Envelope, 
respectively. 
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Fig. 2. — Variation of the CS line profiles for the UCC (left, from top to bottom: Models HCi, HC3 and 
HC5) and the Envelope (right, from top to bottom: Models HC2, HC4 and HCg) when various internal hot 
core structures are used. From top to bottom, we have applied no profiles, only a profile in density, and a 
profile in both density and temperatures (see Sect|3]and SubsectfCTJ). The CS line profiles are expressed 
in antennae temperature (in K) using the IRAM and the JCMT resolutions for low-J and high-J CS lines, 
respectively. 
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Fig. 3. — Variation of the CS line profiles for the UCC with time (age plotted in the third dimension: Models 
HC 7 , HC 8 , HC 5 , HC 9 , HCio, HC12 and HC13, respectively). 
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Fig. 5. — Variation of the CS integrated line fluxes for the UCC (right hand side) and the Envelope (left 
hand side). See text in Subsect. 14.31 
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Table 1: Summary of model runs. "-" means no density or temperature profiles whereas "+" means a profile has been applied. The CS line 
fluxes are displayed in units of 10~ 15 Wm~ 2 . For each line, we specify the corresponding frequency in GHz. We have giving examples of 
paired models (see brackets in columns 1 and 2) representing the coupling of the UCC with the Envelope for a particular 
age. In column 9 we have listed only the values of the CS abundance found at the outer and inner shells for each model. 
The evolution of the CS abundance alongside the core can thus be studied by looking continously at the columns 9 of the 
models belonging to the same pair. 
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I Tin 

HC 8 


UCC 


0.03 


+ 


+ 


O fl \y 

Z.V X 


1U 


O T/" 1 A 00 > , 1 A — 9 

3.76 — 9.82 x f 


r* r* ia — 5 

3.56 x 10 


fi OA 1A — 3 

6.39 x 10 


9.22 x 


10 


1 /IO 1A — 1 

1.42 x 10 


I I Tin 

HC 9 


UCC 


0.03 


+ 


+ 


5.0 x 


fO 3 


0.82 — 3.71 x fO 


1 1 r\ 1 a — 4 

1.10 x 10 


1.52 x 10 


l.Oo X 


IU 


2.42 x 10 


rl^io 


TTPP 


U.UO 


1 

t 


1 


7.0 x 


fO 3 


1.10 — o.oy x iu 


1 7Q y 1 A— 4 


±.yo x iu 


1.73 x 


10- 1 


Z. 1 ^4 X IU 


HC n 


UCC 


0.03 


+ 


+ 


f.O x 


fO 4 


1.64- 9.33 x f0~ 8 


2.61 x 10~ 4 


2.35 x 10" 2 


1.88 x 


10- 1 


3.00 x IO" 1 


HC12 


UCC 


0.03 


+ 


+ 


5.0 x 


fO 4 


0.55 - 4.92 x 10~ 7 


7.59 x 10~ 4 


3.28 x 10~ 2 


2.29 x 


10- 1 


3.64 x IO" 1 


HC13 


UCC 


0.03 


+ 


+ 


f.O x 


fO 5 


8.89 - 4.79 x fO" 7 


9.37 x 10~ 4 


3.48 x 10" 2 


2.39 x 


10- 1 


3.74 x 10" 1 


HC14 


Env. 


0.f5 


+ 


+ 


f.O x 


fO 3 


3.42 - 2.62 x 10~ 9 


2.66 x 10~ 3 


2.79 x 10~ 2 


5.52 x 


10- 2 


5.60 x 10~ 2 


HC15 


Env. 


0.f5 


+ 


+ 


2.0 x 


fO 3 


3.81 - 2.44 x 10~ 9 


2.86 x 10~ 3 


2.83 x 10~ 2 


5.59 x 


10- 2 


5.72 x 10" 2 


HC16 


Env. 


0.15 


+ 


+ 


5.0 x 


fO 3 


5.f5-2.50x 10~ 9 


3.47 x 10~ 3 


2.97 x 10" 2 


5.86 x 


10- 2 


6.17 x 10~ 2 


HC17 


Env. 


0.15 


+ 


+ 


7.0 x 


fO 3 


6.f6-2.83 x 10~ 9 


3.87 x 10~ 3 


3.06 x 10~ 2 


6.04 x 


10- 2 


6.47 x 10~ 2 


HCis 


Env. 


0.15 


+ 


+ 


f.O x 


fO 4 


7.74-3.51 x 10~ 9 


4.42 x 10~ 3 


3.17 x 10" 2 


6.28 x 


10- 2 


6.85 x IO" 2 


HC19 


Env. 


0.15 


+ 


+ 


2.0 x 


fO 4 


13.77- 7.29 x 10~ 9 


5.83 x 10- 3 


3.45 x 10~ 2 


6.88 x 


10- 2 


7.73 x 10~ 2 


HC20 


Env. 


0.15 


+ 


+ 


5.0 x 


fO 4 


2.79 - 2.06 x 10~ 8 


7.45 x 10~ 3 


3.82 x 10~ 2 


7.70 x 


10- 2 


8.78 x IO" 2 


HC21 


Env. 


0.15 


+ 


+ 


7.0 x 


fO 4 


3.45 - 2.58 x 10~ 8 


7.90 x 10~ 3 


3.92 x 10" 2 


7.94 x 


10- 2 


9.05 x 10~ 2 


HC22 


Env. 


0.15 


+ 


+ 


f.O x 


fO 5 


4.17-3.21 x 10~ 8 


8.27 x 10~ 3 


4.03 x 10~ 2 


8.18 x 


10- 2 


9.31 x 10~ 2 


HC23 


Env. 


0.15 


+ 


+ 


5.0 x 


fO 5 


9.53 - 7.29 x 10~ 8 


9.79 x 10~ 3 


4.58 x 10" 2 


9.43 x 


10- 2 


1.06 x 10" 1 


HC24 


Env. 


0.15 


+ 


+ 


1.0 x 


fO 6 


11.89- 1.33 x 10~ 8 


1.00 x 10" 2 


4.61 x 10" 2 


9.48 x 


10- 2 


1.05 x IO" 1 



